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WALL INTERFERENCE IN WIND TUNNELS WITH SLOTTED AND 

POROUS BOUNDARIES AT SUBSONIC SPEEDS 

By Barrett S. Baldwin, Jr., John B. Turner, 
and Earl D. Khechtel 


SUMMARY 


Linearized compressible-flow analysis is applied to the study of 
wind-tunnel-wall interference for subsonic flow in either two- 
dimensional or circular test sections having slotted or porous walls. 
Expressions are developed for evaluating blockage eind lift interference. 


INTRODUCTION 


In solid-wall wind tunnels the effects of blockage severely limit 
model sizes that can be tested at high subsonic . speeds j in fact, the 
model must become vanishingly small as sonic speed is approached. It 
has been demonstrated that if the walls are ventilated (e.g., slotted 
or porous) then blockage is reduced and much larger models can be 
tested. However, wall-interference effects, although reduced, still 
exist and must be evaluated in order to correct the wind-tunnel data to 
free-air conditions. 

It is the objective of the present investigation to analyze two of 
the principal wall- interference effects, blockage and lift interference, 
for two- and three-dimensional subsonic flows in ventilated test 
sections, where blockage refers to the mean incremental velocity induced 
in the vicinity of the model by wall interference and lift interference 
is the mean upwash so induced. In the three-dimensional case it is 
convenient to perform the analysis for a circular test section. The 
results obtained for the circular test section may be applied to a square 
test section of equal cross-sectional area since the wall interference 
at the center of the tunnel should be relatively insensitive to such a 
change in the shape of the wall. 
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SYMBOLS 


factor in Fourier integral transform of 
slot width of slotted wall (see fig. l) 
wing span of model wing 

constant factor in nonlinear term of boundary equation 
Fourier integral transform with respect to x of cp 
dummy variable of Fourier transform 
half tunnel height 

modified Bessel function of the first kind and order zero 

modified Bessel function of the first kind and order one 

modified Bessel fiinction of the second kind and order zero 

modified Bessel function of the second kind and order one 


slot constant, — — In 

' It 


sin ( a 


slot separation of slotted wall (see fig. l) 
lift on the model 
free-stream Mach number 

T 

parameter proportional to size of two-dimensional model, ~ 

parameter proportional to size of three-dimensional model, ~ 

coordinate in the direction of the outward normal to the wall 

dummy variable of integration 

porosity parameter 

cylindrical coordinates 

free-streain velocity 


u,v,w 
u*, v*,w* 

^ rwf 

u,v,w 

Au,Av,Aw 

W 

X 

x,y,z 

a 

P 

r 

« 

9 

9i 

9* 

6<p* 


P 

Te 

Tr 

I 


pertxirbation velocity components in the x, y, z directions, 
respectively 

additional velocity components due to the presence of the 
walls 

additional velocity components having rapid spacewise 
variation near the walls 

additional velocity components at the position of the model 
due to the walls 

complex velocity in the y, z plane 

complex variable equal to z + iy (physical plane) 

Cartesian coordinates 

dummy constant in limiting process 

/l - 

circulation 

complex velocity potential 

total perturbation velocity potential, 9i + 9* 

approximate perturbation potential due to model in free air 

additional perturbation potential due to tunnel walls 

additional wall-interference potential arising from non- 
linear term in boundary equation 

free-stream density 

cross-sectional area of two-dimensional model 
volume of three-dimensional model 

function of X, equal to ^ + It) (treinsformed plane) 
Cartesian coordinates in transformed plane 


4 


NACA RM A53E29 


ANALYSIS 

General Statement of the Problem 


The effect of the tunnel walls on the flow around a model, in -the 
case of ventilated walls, can be calculated using the same basic method 
as that used in reference 1 for the closed -wall case. As in reference 1 
the analysis is based on the linearized equation of subsonic compressi- 
ble flow 




^9 ^ = 

^ Sy2 ^z2 ■ 


( 1 ) 


where cp is the perturbation velocity potential of the flow in the 
tunnel. 

Let q) = % + (p», where is the potential of the flow about the 

model in free air and 9* is the potential of the additional flow due 
to the presence of the walls. 

If 9i is taken to be a known solution of equation (l) which 
approximates the true free -air potential at points far from the model, 
9* can be calculated from the fact that the sum 9i + 9* satisfies a 
known bo\indary condition at the wall. Since the values of 9i at the 
wall only are used, any inaccuracy in the value of 9i near the model 
should not affect the calculation of 9* appreciably. 

The primary objective in this procedure is to estimate the change 
in stream conditions caused by the walls at the position of the model. 

It is assumed that the velocity components derived from 9* are con- 
stants near the model which can be subtracted' from the stream velocity 
to obtain the equivalent free air stream velocity. Thus, 

59 * 

Au = — ^ at x=y=z=0 

ox 

is the blockage correction, and 

5cp* 

Aw = at x=y=z=0 

oz 

is the upwash correction in the three-dimensional case. 


Bo\indary Conditions 


In this section a 
boundary conditions of 
will be developed. 


single expression approximately representing the 
solid, porous, and slotted walls and an open jet 
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Let X be the coordinate in the direction of the free stream 
and n the coordinate in a direction perpendicular toJthe x direction. 
Consider a wall which is perpendicular to the n direction (i.e., 
parallel to the free stream) . If the wall is solid, the condition of 
no flow through the wall can be expressed as 

^ =0 at the wall 

on 

In the case of an open jet there is no pressure drop across the jet 
boundary so that there is zero perturbation pressure at the boundary. 

With a disturbance in the stream this boundary does not remain parallel 
to the free stream. However, for convenience, the condition of zero 
perturbation pressure is imposed at a surface parallel to the free 
stream and coinciding with the jet boundary far upstream of the 
disturbance (see ref. 2). Also, for convenience, this surface can be 
called an open wall and the boundary condition can be expressed as 

— =0 at the wall 

Sx 

In reference 3 an average boundary condition for a porous wall is 
derived. The average velocity normal to the wall is assumed to be 
proportional to the press\xre drop through the wall, a linearized 
approximation to viscous flow through a porous medium, and the pressure 
outside the wall is assumed equal to the free-stream pressure. This 
leads to the bomdary equation 

^ + i ^ = 0 at the wall (2) 

Sx R ^n 

The quantity R is a porosity parameter defined by 

Ap = ^ ^ 

R Bn 

where 

Ap pressure drop through the wall 

P stream density 

U stream velocity 

The quantity PU/R can be determined experimentally by measuring the 
mass flow and pressure drop through a sample of the wall under con- 
ditions corresponding to zero stream velocity. 

Porous walls to which equations (2) and (3) are applicable will ♦ 
henceforth be referred to in this report as ideal porous walls. 
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An approximate "boundary equation for a slotted wall is derived in 
Appendix A. The pressure at the slots is assumed constant and equal to 
the free -stream pressvtre. The resulting xmiform boundary condition is 


where K 


^ + K = 0 at the wall 

hx SxSn 

is related to the slot geometry by 


K = - - 2n 


sin 


ia'\l 

- 


ih) 

(5) 


Slotted walls to which equation (U) is applicable will henceforth be 
referred to in this report as ideal slotted walls . 


Solutions for wall Interference based on equations (2) and (4) can 
be obtained in one calculation by combining them in the form 

— + K i ^ = 0 everywhere at the wall (6) 

Sx ^x^n R 3n 

Thus, wall-interference solutions based on eqviation (6) contain, as 
special cases, those of the closed wall (K -> « or 1/R -5> «)^ ideal 
porous wall (K = O), ideal slotted wall (i/R =0), and open Jet 
(K = 0 and 1/R = 0). Furthermore, equation (6) can be assumed to 
describe a slotted wall having mixed potential and viscous flows in the 
slots. In that case the porosity parameter R can be determined 
experimentally, as it is in the case of a porous wall, by measuring the 
mass flow for a given pressure drop through a sample of the wall. 


If it is found that a nonlinear relationship between pressure drop 
and mass flow exists, it may be necessary to add a term of the form 
f(S<P/Sn) to equation (6). This case is discussed in Appendix B. 


In addition to the foregoing interpretations of equation (6) , an 
interpretation identifying it with slotted walls with tapered slots 
(l and a functions of x) and potential flow in the slots is possible. 
This case is discussed in Appendix C where it is found that instead of 
representing viscous effects, the parameter R is related to the taper 
by 


R dx 


X = 0 


(7) 


where K(x) is the same as in equation (5). 


Blockage in a Two -Dimensional -Flow Tunnel 


Under the assumption of infinitesimal model size, the blockage cor- 
rection will be calculated using equation (6) as the boundary equation. 
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Let X be the coordinate in the free-stream direction and y the 
coordinate in the direction perpendicular to the walls. Let h be the 
half tunnel height so that the walls are at y = -h and at y = +h. 

In these coordinates equation (l) becomes 


and equation (6) 




= 0 


becomes 


± K 


5x3y 


R by J 


y = ±h 


=0 


(8a) 

(8b) 


The - signs on the second and third terms are required because at 
the upper wall, n = +y and at the lower wall, n = -y, n being the 
coordinate in the direction of the outward normal to the wall. 


If 9 is replaced by 9i + 9 , equations (8a) and (8b) yield 


8^9^ 




8x^ 8y^ 


(9a) 


89* + 8^9^ + 1 ^ 

8x 8x8y R 8y /y _ 



+ 

8x8y 


1 89i^ 

^ ' y = +h 


(9b) 


These two equations are sufficient to determine 9* when 9i, the dis- 
txirbance due to the model in free air, is known. 


As in reference 1, the disturbance due to the model at zero angle 
of attack in free air is approximated by a two-dimensional doublet 
which can be expressed as 



-H 


( 10 ) 


The reasoning behind the choice of the doublet is as follows: The 

soxxrce-sink distribution representing a nonlifting model contains the 
same total sink strength as total source strength, so that the distant 
flow field would not resemble that of a single source or sink. The 
center of gravity of the source distribution would lie forward on the 
model compared to the center of the sink distribution which would be aft. 
Hence, there would be a dipole or doublet effect of the model at a 
distance. Other higher-order effects would be much less than that of 
the doublet at large distances. The constant 1 % is related to the 
size of the model with sufficient accuracy for present purposes by 
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me = (11) 

where Tg is the area of the model in the xy plane. 

It is convenient to solve equations (9) hy the Fourier transform 
method. The Fourier transform of 9i is defined by the relations 


Gi(g,y) 

9i(x,y) 



(12a) 

(12b) 


Analogous expressions defining the Fourier transforms of 9 are 


G*(g»y) 



^*(x,y) 


00 


-00 


G*(g,y)e'^®'' dg 


(13a) 

(13b) 


With the substitution of equation (l3b) into (9a) and with aji interchange 
in the order of differentiation and integration there results 


r ^ G*(g,y)e‘^^ dg + F ^ G*(g,y)e"^§^ dg = 0 

J ^00 


Performing the indicated x differentiation and collecting the two 
terms under the common integral yields 

-PVg* + dg . 0 


oo 



This integral will be zero if 

-|32g2G*+ = 0 (ll^a) 


Substituting equations (l2b) and (l3b) into equation (9b) and 
proceeding as before yields 
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( . n* - • + 1 ^ ] f ^ n . ^Gl X 1 SGi \ 

I -ig G + ig K — ± - ) = - ( -ig Gi + ig K ^ ± - ) 

\ Sy R ay/y ^ +jj V ay R ^/y ^ 

(l4b) 

Equations (l^a) and (l4b), involving the transforms Gx and G* of 
<Pi and (f>*, can now be solved in the place of solving equations (9a) and 
(9b). The general solution of equation (l4a) is 

G* = A(g) cosh (Pgy) + B(g) sinh (Pgy) 

Substituting the second term of this expression into equation (l3b) 
yields a term of cp* having an odd y dependence. From symmetry 
considerations it can be seen that (p* should have even y dependence. 
Therefore, B(g) must be zero and, hence, 

G* = A(g) cosh (pgy) (l5) 

The unknown factor A(g) will be determined by substituting equation ( 15 ) 
into equation (l4b), but Gx must be known for this purpose. 

In order to find Gx, equation (lO) is substituted into equa- 
tion (l2a) and there results 


Gi (g,y) = H 


X 2 + 


From reference 4, it is found that 


r 3 

2« J X? + 


dx . - 


2P(y| 


^In reference 4, G(g) is defined as G(g) = / F(f)e^^’^^® df. In 


the tables of F(f) versus G(g), F is usually given as a function 
of p where p = i2itf. If P(p) is the fiinction of p given in the 

table, it follows that F(f) = P(p) = P(i2 nf). With the substitution, 

00 

X = 2nf, it is found that G(g) = -i- P P(ix)e^^® dx. Thus, G(g) in 

'Lo 

the tables is the Foxirier transform, as defined in equations (l2) of 
the function of x which results from replacing p by ix . in P(p) . 
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vhere the symbol | \ denotes absolute value. Differentiating both 
sides of this with respect to g yields 


00 

A f 

2n J 


_ . eigXdx = -i-§-e-|g|P|yi 

2nJ_ :?+Fy^ 2 |g 


so that 


Gi = ^ i e'|s|My| 

(g I 

Substituting equations ( 15 ) and ( 17 ) into (l4b) yields 


( 17 ) 


A(g) 


- ig cosh (Pgy) T ig K 3g sinh (3gy) ± r 3g sinh (3gy) 


y = th 


me g r 

|g| L 


ig e"'®l Tig K 




g|p|y 


i (-Me 


JL \ e'lsl^ly 


J y = ±h 


Making the indicated substitutions (y = ±h) and solving for A(g) yields 
the single complex equation 


A(g) = - 


me 

hn 


(it|| - lgK3 .£)a-|s|e>' 


cosh (3gb) + K3g sinh (3gb) + i ^ sinh (3gb) 


(18) 


Substituting equations ( 15 ) and ( 18 ) into equation (l 3 b) yields 


„« me " (l-[|pl8KB cosh (Sgy) e-lB* 

' kn J ; r ^ 

““ cosh (3gh) + K3g sinh ( 3gh) + i sinh (3gh) 

R _ 


<ig 
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Upon separation into real and imaginary parts and making use of the fact 
that an integral from -« to +» is twice the value of the integral 
from 0 to 00 if the integrand is an even function and zero if it is an 
odd function, this becomes 


cosh (Pgy) cos (gx) 


o cosh (3gh) + Kpg sinh (Pgh) 


- sinh (Pgh) 
R 



cosh (3gy) sin (gx) dg 


With the substitution^ q = pgh^ this becomes 
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The integrals indicated in equation (l9) have not been foxmd in 
closed form. However, if the derivatives of 9* are desired at a 
small mamber of points for a small number of values of K/h and P/R, 
the derivatives can be taken in the integrand and the resulting inte- 
grals evaluated numerically or graphically. 

The quantity of primary interest is Au, the value of 
evaluated at x = y = 0. Performing the indicated differentiation on 
equation (l9) and setting x = y = 0, the expression 


Au = - 


me 


4np2h^ 


w 

/ 




.1 1 

-Iv 



cosh (q) + — q sinh (q) 
h 

+ 

C\J 

sinh q 
R 

2 


'qdq 


( 20 ) 


is obtained for the additional stream velocity at the position of 
the model due to the walls. 


Solid wall.- Letting K — > « or 1/R — > « in equation (20) yields 


m. 


Au 


h^ si 


e-<lq 


dq = 


me 

i2t,2 


2np‘="h‘=^J sinh (q) * 24 p^h‘ 

h ° 


Ideal porous wall .- At K = 0 equation (20) becomes 

2 


Au 


h = 0 


p cosh (q) - sinh q 

2it3 h Q r eosh (q) ^ /_ 


I ) sinh (q) 
K 


(21) 


-.2 


e“^qdq (22) 


Ideal slotted vall «- With 1/R = 0 we have 


i=0 

R 


me 

7^ 


I 


e~^^qdq 




,-2l 


(23) 


Open jet .- For both K = 0 and 1/R = 0 

00 

-21, 




1 = 0 

I- 


jtp^h^ 


r 6 ^^*id 
J 1 + e"' 


qdq _ « “le 

•2q M p2jj2 


Au 


(24) 
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In the general case (eq. 20), the value of Au lies between the 
values for the solid wall and the open jet, and Au can be made arbi- 
trarily small by an appropriate choice of K/h and 3/R. Figure 2 is a 
plot of the values of K/h and 3/R at which Au = 0. These values 
were computed nximerically near the ends of the curve and interpolated 
in the middle. Equations (21), (22), and (24) are in agreement with 
the results of reference 3* Figure 3 is a plot of equation (23), show- 
ing the variation of blockage factor with slot parameter for the two- 

dimensional-flow, ideal, slotted tunnel. 

« 

Since the effect on the blockage correction of letting 3 approach 
zero is the same as letting 1/R approach zero, it can be concluded 
that as the stream Mach number approaches unity, the blockage correction 
factor for an ideal porous wall approaches that of an open jet. ■ Simi- 
larly any effect of viscosity or taper of a slotted wall described by 
equation (6) would be suppressed at near sonic speed so that the block- 
age correction factor would approach that of an ideal slotted wall. 


Blockage in a Circular Tunnel 


Again, the blockage due to the wall interference of a very small 
model will be calculated using equation (6) as the boundary equation. 

Let X be the coordinate in the free -stream direction and r the 
cylindrical coordinate perpendicular to the x direction. Let rQ be 
the tunnel radius so that the wall is at r = rQ. Using these coordi- 
nates, equation (l) becomes, in the case of rotational invariance, 

and equation (6) becomes 


( 


5^9 


^ . K 
Sx SxBr 




^ ^^^r = ro 


= 0 


Then 9 


must satisfy the equations 


r2 ^9* 


+ i ( r j = 

r Sr N Sr / 


(26) 


(27a) 


. K £2-% 1 32-* j =- + (27b) 

\ 6x 6x6r R 6r / r = ro \ 6x 6x3r R hr _ 


(See fig. 1.) 
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Again, from reference 1 the free-air solution is approximated by 




'Pi = 7^ 


(x^ + 

The constant mp is related to the size of the model by 


ID. — U 


(28) 


(29) 


where Tp is the volume of the model. Correction factors which take 
into account the shape of the model and the presence of shock waves are 
discussed in reference 5- 

Substituting, as before, the Fourier integral expressions for cpr^ 
and cpi into equations (2?) yields 




(30a) 


( 


, ^ . Bg* .1 Sg* 

-Kxg — 


I • n XT • SGi . 1 

.Klg— *5 — 


r = r. 


r = r. 


(30b) 


The general solution of the first of equations (30) is 
G* = A(g) lo (pgr) + B(g) Ko (pgr) 

where Iq and Kq are the modified Bessel functions of zero order of the 
first and second kinds, respectively (see ref- 6). Since Kq (Pgr) has 
a singularity at r = 0 which would lead to a singularity in cp* at r = 0, 
B(g) must be zero and, hence. 


G* = A(g) Iq (Pgr) 


(31) 


Substituting equation (28) into the Fovirier integral expression of 
9i yields ^ 

X 


„ mr 1 r X igx ^ 

Gi = — — / s-75 e dx 


(32) 


From reference 4 it is found that 


w 

- f 

2 nJ 


y^^Tp^r^ 


e^S^ dx = - Kq ( Pr I g I ) 
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According to equations (l2), the inverse relation 


00 

/ J OHg|) dg = 


J + 3 ^ 1 ^ 

follows. Taking the derivative with respect to x of both sides of 
this equation yields 


(-ig) i Kq (3r|g|)e“^®^ dg = - 


(x= + 


From this, according to equations (l2). 


e^e^ dx 


2ir J (x^ + 


= ig i Kq (pr|g|) 


SO that 


Gi = ^ ig Ko (3r|g| ) 


( 33 ) 


A(g) 


Substituting equations (31) and (33) into equation (30b) yields 
•ig Iq (Pg^o) -ig KPg Ii (Pgro) + - Pg Ii (Pgro) 1 = 

R J 


m-p 


ig 


-ig Ko (Prolgl) -ig K(-p|g|) Ki (Prolgl) + 


^ (-p|gD Ki (pro|g|) 


the derivatives of I© and being found in reference 6. Upon 
solving for A(g), substituting the result into equation (l3), and 
substituting q for pr^jg there results 


Ki(q)lo(q) + Ko(q)li(q) Uo ( ^ )cos )q dq 



Ideal porous wall.- At K = 0 equation ( 35 ) becomes 


Au 


( 37 ) 


_ “r 2 

00 

r. 

Ko(q)lo(q) - ^ Ki(q)li(q) 

q^dq 

$- = 0 
^0 

> 

io(q) 

2 

+ 

^ Ii(q) 
Lr J 

2 


Ideal slotted vail .- With (i/r) = 0 we have 



Open jet .- For both (i/r) = 0 and K = 0 


Au 


my 


2 rKjq) 


P = 0 

R 

Tr, 




q^dq = 


nir 


0.63 




( 38 ) 


( 39 ) 


The values of Au in equation ( 35) lie between the values for the 
closed wall and the open jet. Figure is a plot of values of K/vq 
versus P/R at which Au =0; the shape of the curve was calculated neajr 
the ends and interpolated in the middle. A graph of equation ( 38 ) appears 
in figure 5 showing the variation of blockage factor with slot parameter 
for the cylindrical, ideal slotted tunnel. 

Again, letting p 0 has the same effect on the blockage cor- 
rection as letting (i/r)->0, so that near sonic speed the ideal porous 
wall should act like an open jet and the slotted wall should act like 
an ideal slotted wall. 


Lift Interference in a Circular Tunnel 


The upwash correction will be calculated using the infinitesimal 
model size approximation and the approximate boundary condition of 
equation (6). 

Let X be the coordinate in the free-stream direction, z the 
coordinate in the direction of lift on the model, and y the remaining 
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Cartesian coordinate. The cylindrical coordinates r and B are related 
to the Cartesian coordinates hy r + and z = r sin 0. It can 


be seen from reference 1 that the appropriate free air solution is 


m Tb 


1 + 


sin_0 




(40) 


which is the potential of a horseshoe vortex having infinitesimal span. 
The fact that the actual span is finite introduces higher-order terms 
which are negligible at distances large compared to the size of the 
model. Here Fb is related to the lift on the model by 


L = pUFb 


(4l) 


The Fourier transform with respect to x of 9i cannot be found. 

An arbitrary parameter a will be introduced into the potential so that 
the Fourier transform of a related function can be found. Let this 
fxmction be 


cPi- = 


Fb 

4jt 






S -CLt/ x^ + 3^r^ 


a 


8x 


sin 0 (L 2 ) 

r 


so that 


lim 9x' = 
a — > 0 


Then a will be eliminated from the resulting 9*"* by taking the limit 
with a at zero. 

From reference 4 it is found that the Fourier transform of 


, g g|3r Ki * c?) 


+ or 


. By the use of an inversion. 


differentiation, and reinversion, as in the derivation of equation (33)> 


it is found that the transform of + 3^^ 

dx 


IS 
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lga3rKi(3ryg^+ g^ ) 


so that the transform of 9i* is 


Gi' = M (a+ig) Ki(3ryg^-H3.^) sin 0 

/y/ g2 + a2 


(43) 


From physical reasoning, it can be seen that the sin 9 factor 
of cpi will result in the same angle dependence for cp*, so that sin 6 
will be a factor common to cpi', cp*' and their Fourier transforms Gi' 
and G* . With cylindrical coordinates, equation (l) is 


p2 ^ + 1 A/'r Y 1 5^9 _ 

r 3r \ J 


In the case of sin 6 dependence, this becomes 


+ lA(^r^).i_cp = 0 

r Sr \ Sr / r^ 


and the boundary condition remains the same as in equation (26) . Sub-, 
stituting the Fourier integral expressions for cpi^ and cp*’ in these 
equations, as before, yields 


-pVg*’ 


r ^ J 



.(41^a) 


f-lgG*'- Kig Al 

V Sr 


+ 1 

R Sr ' psrp 



-Kig 

Sr 


1 SGi» 
R Sr 


r=r 


o 


(44b) 

The solution of equation (44a) which has no singularity at r = 0 is 


G*’ = A(g) Ii(prg) sin 0 


(45) 



Replacing Ii (qr/ro) by its pover series expansion, differentiating with respect to 
which is equal to r sin 9, and setting r = 0 yields 
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Ideal porous wall » - At K = 0 equation (48) becomes 


Awl 


Tb 


K 


4jtro^ 


•=o 


1- P i 

R It 


[Ki(q)lo(q) + Ko(q)li(q) ] q^dq 




[Ii(q) - qlo(q)]' 


Ideal slotted wall .- When 1/R is set equal to zero in equa- 
tion (48) y a limiting process is required to obtain the correct result 
at (p/R) = 0. The result of this process is 



Open jet >- Letting K = 0 in equation ( 50 ) yields 


(50) 


Av 


3 

r 


= o 


K 


^o 


o 


Fb 

4nr/ 


The value of Aw in the general case of equation (48) lies between 
the values for the closed wall and the open Jet. 


CONCLUDING REMARKS 


A method of evaluating wall interference of partly open walls involv- 
ing mixed potential and viscous flows has been presented. Expressions for 
blockage and lift interference for both slotted and porous walls have been 
derived. Some new details of the method may prove useful in other theo- 
retical treatments of this type of problem. 

The results of the analysis indicate that near sonic speed, the 
blockage correction for an ideal porous wall approaches that of an open 
Jet. Similarly, any linear viscous or taper effect of a slotted wall is 
suppressed near sonic speed, so that the blockage correction approaches 
that of an ideal slotted wall. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., May 29, 1953 
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APPENDIX A 

DERIVATION OF THE BOUNDARY EQUATION 


In this appendix an approximate smoothed or average boundary equa- 
tion for a slotted wall will he derived. 


An ideal slotted wall has zero perturbation pressure at the slots 
and zero normal flow at the strips. These conditions can be expressed 
as 


= 0 at the slots 

Sx 


^ = 0 
Sn 


at the strips 


(Al) 


When the slot spacing and model dimensions are small compared to 
tunnel dimensions, the perturbation flow can be separated into a rapidly- 
varying and a relatively uniform part so that the two parts can be 
investigated separately. It will be shown that the effect of the rapidly 
varying part can be replaced by a condition on the relatively uniform 
part. 


Let cp, u, V, and w represent the rapidly varying part of the flow 
field and <p, u, v, and w the remaining part of the perturbation flow. 
For a plane -wall at z = h equations (A1) require that 


u 

+ u = 0 

at the slots 

w 

+ w = 0 

at the strips 

In addition 






U = V 

= w = 0 

far from the wall 


(A2) 


(A3) 


In order to solve for cp, use can be made of the fact that u, v, 
and w are nearly constant at the wall compared to u, v, and w, so 
that u and w can be considered constant in equations (A2). 


Since the slots lie along the x direction, u is nearly constant 
in the x direction, so that 3{3ySx can be neglected compared to 

Sv/Sy and ^/Sz. As in slender airplane theory, this leads to a two- 
dimensional crossflow for which 


2k 
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(A4) 


Then $ can be of the form Q(x) f(y, z) where Q(x) is a slowly- 
varying function of x, and f(y, z) satisfies the two-dimensional 
Laplace equation . Then, 


u = ^ ^ f(y,z) 

ox ox 

and since u is equal to -u, a constant, at the slots, f(y, z) must be 
constant at the slots or 

V = = Q ^ = 0 at the slots 

dy oy 


This equation can be used to replace the first of equations (A2), 
and altogether there results 

V = 0 at the slots 

w = -w = constant at the strips f (A5) 

^ = V = w = 0 far from the wall 

assuming that BQ/Sx is not absolutely zero. 

Equations (aU) and (A5) can be solved using the conformal trans- 
formation technique. 


Let the wall at z = h be slotted periodically from y = -« to 
y =. +00 with the center of a slot at y = 0. Let a be the slot width 
and Z the slot separation (see fig. l) . It is sufficient to consider 
only one period of the periodic flow configuration so that solid bounda- 
ries can be placed at y = ±(z/2) and attention confined to the region 
between them. 

Let X = z + iy be the complex physical plane and I = C + the 
transformed plane. Let be the complex velocity potential such that 
W = w - iv, the complex velocity, is equal to dO/dX, and 9 = real part $. 
Then equation (a 4) is satisfied if is any analytic function of X. 

The boundary conditions (A5) are satisfied by finding the analytic fionc- 
tion I (X) which transforms the boundaries in the X plane into a con- 
figuration for which the flow field with the desired flow at the bound- 
aries can be found. 
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In the X plane the region under consideration is from z = -«* 
to z = h and from y = -(z/2) to y = +(l/2). The slot lies on the 
line X = h + iy from y = - (a/2) to +(a/2) . • (See fig. 6.) The 
two half strips lie on the same line from y = -(z/2) to y = -(a/2) 
and from y = (a/2) to y = (l/2). The remainder of the solid boundary 
lies along the line X = z — 1 (z/2) from z = -»to z=h and along 
the line X = z + i(z/2) from z = -» to z = h. 

The transformation which will place the origin at the wall, and the 
domain londer consideration in the right half -plane is 


li = h - X (a6) 

The effect of the transformation on the positions of the boundaries is 
obtained by substituting the equations representing the boundaries into 
equation (a 6) • Thus, the line X = h + iy becomes + i^i = -iy 

in the plane and the slot lies on this line from t|x = -(a/2) to 

Tlx = + (a/2). The two half strips lie on the same line from t)i = -(z/2) 
to -(a/2) and from a/2 to Z/2. Similarly, the remainder of the solid 
boTindary in the plane lies on the line = h - z - i(z/2) from 

^ 1=0 to and on the line = h - z + i(z/2) from Ci = 0 

to ^ 1 = 00 . 


To satisfy the boundary condition at the strip, a term w|^ = w(h-X) 
is added to the potential in the transformed plane. 

The remaining boundary -value problem is solved with the aid of two 
further transformations. 


The transformation 



(a7) 


transforms the region under consideration in the li plane to the entire 
right half of the I 2 plane, as can be seen by following the procedure 
outlined for the first transformation. It is found that in the I 2 
plane, the slot lies along the imaginary axis from -i sin (jta/2Z) to 
+i sin («a/2Z) and the solid bovindary including the two half strips lies 
along the imaginary axis outside of ±i sin (ita/2Z). 

The transformation 

= iz sin® 




(A8) 
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transforms the region under consideration to the positive real half of 
the i plane, excluding the circle of radius sin (na/2Z) centered at 
the origin. The slot lies along the half of this circle in the right 
half plane, and the complete solid boundary lies along the imaginary 
axis outside the circle. Thus, a source at the origin in the | plane 
will satisfy the condition of no flow through the solid boundary and no 
flow across the slot. Hence, the desired potential is 

$ = AZn ( I) + wli (A9) 

Substituting equations (a 6) , (AT), and (a 8) in (A9) yields 


= AZn -Isinh 




I (h-X) 




sinh^ 


1 (h-X) 



+ w(h-X) 


The constant A is evaluated by the last of equations (A5), 
being the equivalent of that equation. 

As X -®, the potential simplifies to 


=0 

X 


<D 


= AZn |2 sinh | | + ^ (h-X) 

t f (h-X)' 


= AZn 


+ w(h-X) 


= A i (h-X) + w (h-X) 

L 

The neglected terms are of order smaller than l/h-X. Thus, will be 
zero at infinity if A = -(z/jt)w so that 


® = - w 


Zn -jsinh ^ (h-X) + /^inh^ ff (h-X) 


+sin 


(AlO) 


Then 


ly=o 


= - w ^ Zn 


sinh 




+ y sinh^ (h-z)J + sin^^^'^1, + w(h-z) 
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It was assumed that the velocity potential is constant over the slot. 
The value of this constant is 


z=h 


0 ' ■“ « (i) . 


= 9 


at the slot 


and 


sin 


«a 

21 


= u at the slot 


b<p _ ^ r Z 

Sx Sx \ ^ 

Substituting this in the first of equations (A2) yields 


u + 


sin 


1 ^ _ 
J J Sx " 


0 at the slot 


Since it was assiomed that u and w do not vary appreciably from 
slot to strip, this equation applies everywhere at the wall and yields 


for a plane wall. 



at the wall 


(All) 


In considering a curved cylindrical wall, it appears that the above 
results are not altered appreciably if the radius of curvature of the 
wall is everywhere large compared to the slot spacing.^ Hence, it can be 
assumed that equations (All) are applicable to any slotted wall. 


^For a circular cylindrical slotted wall, solution of the boundary value 
problem 

a? 

— = 0 at slots 

he 

— = -w at strips 

Sr 

' 9=0 at r = 0 

yields a value for K identical with that obtained from equation (All) . 
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In an attempt to take viscous effects in the slots into account, it 
can be assumed that as in the case of the "ideal" porous wall there is a 
pressure drop through the wall which is proportional to the normal 
velocity at the wall. In that case the first of equations (A2) is 
replaced by (u + u) +(i/R)w = 0 at the slots which leads to 

+ K = 0 everywhere at the wall (A12) 

8x oxon R on 

where K remains the same as in equation (All) and R is to be deter- 
mined experimentally. 

Equation (AlO) can be used to calculate the neglected variations 
in flow quantities near the wall and also at the model if the variations 
are not negligible there. 


DISCUSSION OF A CRITICAL ASSUMPTION IN THE DERIVATION 
OF THE BOUNDARY EQUATION 


It is assumed that the perturbation pressure at the wall is pro- 
portional to u + u. This is a good approximation only if 

(v + + (w + w)^«2 (u + u) U 

in addition to the usual requirements for linearization. This addi- 
tional requirement can be reasonably relaxed to 

(w + w)^ «2uU at the wall (A13) 


Equation (AlO) indicates singularities in w at the edges of the 
slots. Experience with wing leading edges indicates that this dis- 
crepancy can be reasonably Ignored. However, equation (A13) should at 
least be satisfied at the center of the slots. The effect of this 
condition on u and v will now be determined. Differentiating 
eqiiation (AlO) yields 


cosh 


W = 


d® 


rt 

I 


= w 


/ 


sinh'" 


7 (h-x) 

6 


(h-x) 


+ sin 


— w 


/ Tta \ 

1 , 2 !) 



4 . 
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Then 



so that equation (A13). "becoines 




« 


sin^ 



2Uu 


or 


« 2 sin^ 



at the wall 


(Al4) 


This result places a lower limit on the ratio of open to total area 
(a/l) for which the results of this analysis can he expected to apply 
to slotted sections. 
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APPENDIX B 

BLOCKAGE IN A CIRCULAR TUNNEL WITH 
NONLINEAR VISCOUS EFFECT 


If, in the e 3 qperimental determination of the porosity parameter R, 
it is found that an additional term of the form f (Scp/Sn) is needed in 
the boundary equation, the result is 

^ + K i — + f f"— ^ = 0 at the wall (b1) 

dx^n R Bn \^n/ 

or 

^ + + = (B2) 

Bx BxBn R Bn Bx BxBn R Bn \ Bn Bn/ 


This type of equation cannot be solved exactly by any presently known 
method because of the nonlinearity. However, the equation becomes linear 
if the <p* in the nonlinear term is neglected under the assumption that 

f _ f < < f at the wall (B3) 

\ Bn Bn / \ Bn / Bn y 


Hence, the equation 

atthewall 

Bx BxBn R Bn Bx BxBn R Bn \ Bn/ 

(Bk) 

will be discussed. 

The transforms G* and Gi of 9* and 9i remain the same as before, 
with the exception that the factor A(g) must be evaluated from the 
transform of equation (b 4) . For this purpose, the transform of 
f(B9i/Br) is needed and can be found by evaluating 

e^®^ dx (B5) 



" 00 


r=ro 
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The transform of equation (b 4) is found to be 


A(g) - iglo (pgr^) - igKpgli + i pgli 
L R 


- ^ ig (^-igKo + ig P [g [KKi - J 3 Is “ Gf 


Upon solving for A(g) and substituting in the inverse transform 
expression for T* equation (3^) is obtained plus the additional term 


oo 


J Pgli(Prog) - ig [lo(Prog) + KPgli(prog) 


GfIo(Prg)e~^SX^g 


With the substitution q = Prog this becomes 

1 


Scp* = - 


I iM- 


i flo(q) + ^ qli(q) 

To 


• 

> Gfe" O-o M (b6) 


As an example of the use of equation (b 6) let 

3 


f { ^ 

\3n/ \5n/ 


(B7) 


This fvinction is chosen rather than one proportional to (S9/Sn)^ 
because the pressure drop should be an odd function of the normal 
velocity. 


Hence, 



OP^ro)"" 



X 


3 


[ x^ + P^ro^] 


15/2 


eig^dx 


From reference 4 it is found that 
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r 1 ^igx ^ ^ IgT K 7 (pro|g|) 

2« J. ' 'j^v (2Br„)^ 

Differentiating three times with respect to g yields 



Substituting this in equation (b 6) and separating real and 
imaginary parts yields 
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Upon differentiation with respect to x and setting x=r=0, there 
results 


8Au = - 3 ® — 


mr 


Y 3^ 3 r 

\Un33ro®/ I 


d® 

<l‘^K7(q) 

dq 

^ Il(«l) 

LB J 

2 

+ 

loCq) + Ii(q)1 

^o 


From reference 6 it is found that 


and 


q‘^K7(q) = (1152 + l68q^ - 5q^)q^Ki(q) + (576 + 12q^ - q‘‘)q®Ko(q) 


3® 12 



SO that 


P 12^ 

1 r 

(1152 + i68q^-: 

5 q^)q^Ii(q)Ki(q) dq 

R 5005 n ’ 

j 

0 

^ Ii(q) 

2 

+ 

io(q) + — ii(q)l 

ro 1 



/ 


(576 + 12 q 

- q^) q®Ii(q)K^(q) 

j !.(.)■ 

2 

+ 

’lo(q) + zr Ii (qr 
^ 0 

2 


dq 


(BIO) 


Setting K/tq equal to zero in (BIO) and recombining it with 
equation (37) yields 


(q)Ko(<l) - Ki(q)li(q) 
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compared with the value of 
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Since the factor r-^/knvQ^ is very simll, the nonlinear contribu- 
tion to the blockage correction should be negligible at subsonic speeds 
where 0 is not small. But as contrasted to the linear viscous effect 
which becomes small as M approaches 1^ this nonlinear viscous effect, 
may become large. 
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APPENDIX C 

DERIVATION OF THE BOUNDARY EQUATION 
FOR A TAPERED SLOTTED WALL 


In this appendix an approximate smoothed boundary equation for a 
slotted wall with tapered slots will be derived. 


The development in this case is identical with that in Appendix A 
up to the point where it is found that 



at the slots 


or 

9 = w K at the slot 


(Cl) 


When the slots are tapered, the slot parameter K is a slowly 
varying function of x which can be expanded in a power series about 
X = 0* so that neglecting the higher order terms 


X 

X = 0 

Then 


K(x) = K 


X = 0 


dx 


(p = K 


w + 


X = 0 


dx 


xw at the slots 


X = 0 


and 


5x 


bw ^ 6K 


X = 0 Sx dx 


dK 
w + — 

X = 0 ^ 


X = 0 


bw 

bx 


X at the slot 


Assiiming that the last term will be negligible compared to the others, 
we have 


u = K 


X = 0 


bw dK 
^ dx 


w at the slots 


(C2) 


X = 0 
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Substituting this equation in equation (A2) of Appendix A yields 


u + K 


X 


0 Sx dx X = 0 


at the slots 


Then, as in Appendix A, since u and w do not vary appreciably from 
slot to strip 


Sx 


X = 0 ^xSn dx 


X = 0 


^=0 


everywhere at the wall 


K(x) 


^(x) 

It 


Zn 


sin 


na(x) 


} (C3) 


j 


Eence^ the parameter 1/E in equation (6) of the text can be interpreted 
as the quantity dK/dx | ^ _ q '^hich is related to the taper of the slots 

in a slotted vail having purely potential flow in the slots. 


38 


NACA RM A53E29 


REFERENCES 


1. Baranoff, A. v.; Tunnel Correction for Compressible Subsonic Flow. 

NACA TO 1162 , 1947 . 

2. Glauert, H.: Wind Tunnel Interference on Wings, Bodies and Air- 

screws. R&MNo. 1566 , British A.R.C., 1933« 

3. Goodman, Theodore R.: The Porous Wall Wind Tunnel. Part II - 

Interference Effect on a Cylindrical Body in a Two-Dimensional 
Tunnel at Subsonic Speed. Cornell Aeronautical Laboratory, Inc., 
Rep. No. AD-594-A-3, 1950. 

k. Campbell, George A., and Foster, Ronald M.: Fourier Integrals for 

Practical Applications. D. Van Nostrand Company, New York, 1948. 

5. Gothert, B.: Wind-Tunnel Corrections at High Subsonic Speeds 

Particularly For an Enclosed Circular Tunnel. NACA TO 1300, 1952. 

6. McLachlan, N.W.: Bessel Functions for Engineers, Oxford University 

Press, London, 1946. 


NA.CA FM A53E29 


39 



(a) Two -dimensional -flow slotted test section. 



Figure /. - Cross -sectional diagrams of the two slotted test sections in 
a plane perpen dicula r to the free stream. 
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Figure 2. - Simultaneous values of sJot parameter and porosity 
parameter for zero blockage in o two-dimensionai-fiow tunnel. 



Figure 3r Blockage factor as a function of slot parameter in a 
two-dimensional-flow, ideal, slotted tunnel. 
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Figure 4.- Siimtltaneaus values of slot parameter and porosity 
parameter for zero blockage in a circular cylindrical tunnel. 



Figure 5 .* Blockage factor as a function of slot parameter in 
d^ircular cylindrical» ideals slatted tunnel. 



Transformed p!ane,(^ 


Transformed plane, S 


Figure 6. 


- Sketch of complex planes (Appendix A). 
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